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Abstract 

The breaking down of the equivalence principle, when discussed in the con- 
text of Sikorski's differential space theory, leads to the definition of the so-called 
differentially singular boundary (d-boundary) and to the concept of differential 
space with singularity associated with a given space-time differential manifold. 
This enables us to define the time orientability, the beginning of the cosmo- 
logical time and the smooth evolution for the flat FRW world model with the 
initial singularity. The simplest smoothly evolved models are studied. It is 
shown, that the cosmological matter causing such an evolution can be of three 
different types. One of them is the fluid with dark energy properties, the sec- 
ond the fluid with attraction properties, and the third a mixture of the other 
two. Among all investigated smoothly evolved solutions, models qualitatively 
consistent with the observational data of type la supernovae have been found. 

PACS: 04.20, 04.60, 02.40 

Key words: singularities, dark energy 

1 Introduction 

It is generally believed that the Universe had a beginning and everything indicates 
that it indeed had the beginning. According to contemporary ideas the Planck era 
is the beginning of the Universe. During the flow of cosmic time the known Universe 
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emerges from the Planck era. Certainly it is not an immediate process and one can 
imagine that individual space-time structures emerge from the Planck era gradually. 
One can expect that the simplest structures such as a set structure or a topology on 
this set can emerge first. The appearance of the manifold structure is the key event 
of the process. Then the equivalence principle, in its non metric form, appears and 
makes possible a creation of " higher" structures such as, for example, the Lorentzian 
structure. From this perspective it is interesting to ask how the space-time geometry 
is " nested" in theories more general than the differential geometry since the emerging 
process of gravitation has to be associated with a theory which is more general than 
the geometry of space-time manifold. 

In this paper the Sikorski's differential spaces theory (see Appendix A) is applied 
to the discussion of the breaking down process of the equivalence principle. Within 
this theory the non-metric version of the equivalence principle is one of the axioms 
of the space-time differential manifold (d-manifold) definition. Let us consider points 
at which this axiom is not satisfied. The set of this points forms the so-called differ- 
entially singular boundary (d-boundary). This issue is discussed in Section 2. The 
construction of this type of boundary is also described. 

In Section 3, the construction from Section 2 is used for building the differential 
space (d-space) with d-boundary for the flat FRW model. The issue of the prolon- 
gation of cosmic time and time-orientability to the d-boundary for this model is dis- 
cussed in Section 4. In Section 5, we describe the concept of the smooth evolution of 
this model starting from the initial singularity. We find the simplest smoothly evolved 
flat FRW models in Section 6. Some of them (see Section [7]) have the kinematical 
properties which are qualitatively agreeing with the recent cosmological observations 
Pfl|2]. Namely, these models begin their evolution with the decelerate expansion which 
changes into the accelerated one. 

The discussion of properties of the cosmological fluid which causes that type of 
smooth evolution for the simplest cosmological models is carried out in Section 7. 
Behaviour of any smoothly evolving flat model in a neighbourhood of the initial 
singularity is investigated in Section 8. Finally, in Section 9, we summarize the main 
results of this paper 

2 Sikorski's differential spaces and singularities 

Space-time is a 4-dimensional, Lorentzian and time-orientable d-manifold of the class 
C°° |3J. This definition is a mathematical synthesis of what is known about gravita- 
tional field. It implicitly includes the equivalence principle [U EJ E] • 
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The equivalence principle, in its non-metric version, is implicite in the axiom 
stating that space-time is a d-manifold which means that it is locally homeomorphic 
to an open set in R n , where n — 4. Differential properties of these homeomorphisms 
(maps) are determined by the atlas axioms [7| through the assumption, that the 
composition of two maps (pi o cp^ 1 , as a mapping between open subsets of R n , is a 
diffeomorphism of the class C k , k G N. In the present paper it is assumed that these 
diffeomorphisms are smooth. It is worth to notice, that the non-metric version of the 
equivalence principle is encoded into two levels of the classical d-manifold definition, 
namely in the assumption of the existence of local homeomorphisms to R n and in the 
axioms of atlases. 

One can look at space-time, or more generally at a d-manifold, from the viewpoint 
of theories that are more general than classical differential geometry. In these theories 
d-manifolds are special examples of more general objects. Examples of theories of this 
type are: the theory of sub-cartesian spaces by Aronszajn and Marshall [H |5j, the 
theory of Mostow's spaces [10] or Chen's theory [II] . In the present paper we study 
the initial singularity problem and the problem of the beginning of cosmological time 
with help of Sikorski's differential spaces (d-spaces for brevity) [TJl [321 [33]. In this 
theory, the generalization level does not lead to an excessive abstraction and therefore 
the d-space's theory may have applications in physics. 

Generally speaking, d-space is a pair (M, Ai), where M is any set, and Ai is a 
family of functions ip: M — > R (see Appendix EJ. This family satisfies the following 
conditions: a) Ai is closed with respect to superposition with smooth functions from 
C°°(R ri , R), and b) Ai is closed with respect to localization. The precise sense of these 
notions is not important at the moment. The d-space (M, Ai) is also a topological 
space (M, 7>f), where tm denotes the induced topology on M given by the family 
Ai. The topology Tm is the weakest topology in which all functions from Ai are 
continuous. With in the theory of d-spaces the d-manifold definition assumes the 
form 

Definition 2.1 Non empty d-space (M,Ai) is n- dimensional d-manifold, if the fol- 
lowing condition is satisfied 

(*) for every p G M, there are neighbourhoods U p G and O G r^n, such that 
d-subspaces (U p ,Aiu p ) and (0,£q^) of d-spaces (M,At) and (R n ,£^) respec- 
tively are diffeomorphic, where := C°°(R n ,R). 

Briefly speaking, d-space (M, Ai) is a n-dimensional d-manifold if it is locally diffeo- 
morphic to open subspaces of R n which are treated as d-spaces (see Definitions IA.8I 
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and IA.9[) . Diffeomorphisms appearing in condition (*) are generelized maps. Unlike 
in the classical definition, these maps are automatically diffeomorphisms. Definition 
12.11 is equivalent to the classical d-manifold definition [T2] . 

In the d-spaces theory, the non-metric version of the equivalence principle is "lo- 
calized" in the single axiom (*). This enables us to give the following definition of 
space-time 

Definition 2.2 A space-time is a non empty d-space (M,Ai) such that 

a) (M,A4) satisfies the non-metric equivalence principle as expressed by the con- 
dition (-k) for n=4, 

b) a Lorentzian metric form g is defined on M , 

c) Lorentzian d-manifold (M,A4,g) is time orientable. 

If the non-metric equivalence principle in the classical definition of space-time is 
thrown out, the time-orientability, the Lorentzian metric structure, the classical dif- 
ferential structure and the topological manifold structure are automatically destroyed. 
Only the topological space structure survives. 

In the case space-time Definition 12.21 the situation is different. Throwing out the 
non-metric equivalence principle in the form of axiom (*) leaves a rich and workable 
structure on which one can define generalized counterparts of classical notions such 
as orientability, Riemannian and pseudo-Riemannian metrics, tensors, etc [121 E31 EH 

nana mi Da esi 123. 

Looking at space-time manifolds from the d-spaces perspective, one can imagine 
a situation when a background of our considerations is a sufficiently "broad" d-space 
(M,A4) which is not a 4-dimensional d-manifold. Let us additionally suppose that 
in (M, Ai) there exists a subset B C M having a structure of a 4-dimensional d- 
manifold satisfying the condition (*) for n = 4, regarded as a d-subspace of (M,JA). 
Naturally, at each point p G B the non-metric equivalence principle is satisfied. 
The points not belonging to B can considered as singular points because at these 
points the four- dimensional non-metric equivalence principle is violated. Among such 
singular points in M the most interesting are accumulation points of the set B in 
the topological space (M, t m ). One can call a set of all singular accumulation points 
of B the singular boundary of (B, Ai b) space-time d-manifold. Singular points in 
M not belonging to the singular boundary can be excluded from considerations as 
unattainable from the interior B of the space-time d-manifold. 

Definition 2.3 Let (M, Ai) and (Mo, .Mo) be Sikorski's d-spaces, such that 
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a) (Mo, Mo) is a n- dimensional d-manifold, 

b) (M,M) is not a n- dimensional d-manifold, 

c) Mo C M and Mq is a dense set in the topological space (M,tm), 

d) (Mo,A4 ) is a d-subspace of the d-space (M, M). 

The set d^Mo := M — M is said to be differentially- singular boundary (d-boundary 
for brevity) of d-manifold (M ,Mq), if for every p 6 d^Mo and for every neigh- 
bourhood Up G r M , the d-subspace (U p ,Mu p ) is not diffeomorphic to any d-subspace 
(V,C°°(R n ) v ) of the d-space (R n , C°°(R n )), where V e r K n. Then, the d-space 
(M,M) is said to be d-space with differentially singular boundary (or d-space with 
d-boundary for brevity) associated with the d-manifold (Mo, Ado)- 

In the following we shall restrict our considerations to the case of Definition 12.31 
for n — 4. The d-boundary definition for a space-time's d-space (M ,Mo) depends 
on the choice of the d-space (M,M). Therefore, a reasonable method of the d-space 
(M,M) construction is necessary. Fortunately, the idea described above suggests 
such a construction. Namely, one chooses a sufficiently "broad" d-space "well" sur- 
rounding the whole investigated space-time M and then one carries out the process of 
determination of all accumulation points for M . These points form the d-boundary 
ddM of space-time. Next, one defines the set M as M := M Uc^Mq. By treating M 
as a d-subspace of the sufficiently "broad" d-space, the d-space with the d-boundary 
(M, M) is uniquely determined. On can easily check that the pair of d-spaces (M, M.) 
and (Mq, M.q), where (M, M.) is determined by above described method, satisfies Def- 
inition 12.31 

The choice of a sufficiently "broad" d-space "well" surrounded the whole investi- 
gated space-time is almost obvious. As it is well known [21], every space-time can be 
globally and isometrically embedded in a sufficiently dimensional pseudo-Euclidean 
space E p ' g , where p, q depends on space-time model. The space E p ' q is also the d- 
space (IR^ 13 , C°°(W +q , K)) which "well" surrounds the studied space-time d-manifold 
(Mo, Mo). ' 

Such a choice of the surrounding d-space is well motivated since all causal prop- 
erties of the studied space-time are taken into account, even if we do not refer 
to them directly. This is because of the isometricity of the embedding. Further- 
more, the concrete form of the isometrical embedding provides us with a practical 
method of constructing the d-structure M for the the d-space (M,M). Namely, 
M = C 00 (M p+( ',M) Af . 
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To test this method of constructing (M, Ai) we shall apply it to the flat FRW 
cosmological model. 

3 A differential space for the flat FRW model with 
singularity 

Let us consider the flat FRW model with the metric 

g = -dt 2 + a 2 (t)(dx 2 + dy 2 + dz 2 ), (1) 

where (t,x,y,z) G W° := D® x M. 3 and D® is a domain of the scale factor . The 
scale factor a : D° 9 t — > ao(t) G M + , a G C°°(D°), is an even real function of the 
cosmological time t. A domain of the scale factor D® C M + , is an open and connected 
set. For convenience, let us assume that t = is an accumulation point of the 
set D® and ^ D a a . In addition, the model has an initial singularity at t — i.e. 
lim t ^ + o,o(t) = 0. In the next parts of the paper, the symbol a is reserved for the 
following function: a: D a — )• E, a{t) = a (t) for t G D®, a(t) = for t = where 

IK-. /;,%.{()}. 

A construction of a d-space with d-boundary for the flat FRW model pass in 
analogical way to the construction of a differential space for the string-generated 
space time with a conical singularity described in [22] . 

Every 4-dimensional Lorentzian manifold can be isometrically embedded in a 
pseudo-Euclidean space |2T] . In particular, the manifold (W°,g) of model ([1]) can 
be isometrically embedded in (M 5 ,-?/ 5 )) by means of the following mapping 

F°: W° — > F°(W°) C M 5 

*?(f , x, y, z) = ±a (t)(x 2 + y 2 + 2 2 + 1) + \ £ 

F°(t, x, y, *) = f a (t)(* 2 + + z 2 - 1) + \ £ (2) 
^3 (*> x ' 2/' z ) = a o(t)x, F 4 °(t, x, y, z) = a (t)y, F 5 °(t, x, y, z) = a (t)z, 

where 77® = diag(—l, 1, 1, 1, 1). 

As every manifold, space-time (W°, g) is also a differential space (W°, W°), where 
the differential structure W° := is a family of local £ ^-functions on W° C M 4 , 
where £ (4) = C°°(M 4 ). On the other hand, the set F°(W°) can be equipped, in a 
natural way, with a differential structure treated as a differential subspace of the 
(M 5 ,£( 5 )), where 8^ = C°°(IR 5 ). Then the family £% w0) of local £ ^-functions 
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on F°(W°) is a differential structure and the pair (F°(W°), 8 F L w0 C) is a differential 
space. In addition, if the integral J Q * dr/do(r) is convergent for every t E D , then 
the mapping F° is a diffeomorphism of the differential spaces (F°(W ),£^, w0 ^) and 
(1V°,>V ) [231 [22]. 

The process of attaching the initial singularity depends on the completion of the 
set F°(W°) by means of points from the surrounding space R 5 in the way controlled 
by the isometry F°. It enables us to define the d-space with d-boundary (with the 
initial singularity) for the flat FRW model. 

Let the mapping F: W — > R 5 denotes a prolongation of F° to the set W := 
D a x R 3 . The values of F are given by formulae (T2J) changing the symbol F° onto 
F . Then for every x,y,z G R the value of F at the initial moment t — is 
F(0, x,y, z) = (0,0,0,0,0), since a(0) = by assumption. Therefore, the cosmo- 
logical initial singularity (d-boundary) distinguished by the embedding procedure is 
represented by the single point: d d F°(W°) = {(0, 0, 0, 0, 0)} C R 5 . 

The pseudo-Euclidean space (R 5 ,?/ 5 )) is a differential space (R 5 ,£( 5 )). The differ- 
ential structure 8^ of this space is finitely generated. The projections on the axes of 
the Cartesian system 7Tj: R 5 — > R, tt^zi, z 2 , Z5) = z i: i = 1, 2, 5, are generators 
of this structure and therefore: 8^ = Gen(7Ti, 7T2, ...,^5). As is well known [151 122] . 
every subset A of a support M of a differential space (M, C) is a differential subspace 
with the following differential structure: Gen(C|^). For finitely generated differential 
spaces (M,C), every differential subspace with a support AcMis also finitely gen- 
erated. Generators of the differential structure on A are generators of the differential 
structure C restricted to A. 

The differential space (F(W),£p) W s) is a differential subspace of the differential 

space (R 5 ,£( 5 )) and represents the d-space with d-boundary for model ([1]). Since 
8^ is finitely generated, the d-structure £p) W y induced on F(W), is also finitely 
generated and 

£f(W) = Gen(7ri|F(w), ^2\f(w), K5\f(w))- 

The d-space (F(W), £p) w \) is not convenient for further discussion because it 
depends on the embedding procedure. Let us define a d-space (W, VV) diffeomorphic 
to (F(W), £p%y\) which will enable us to apply the Sikorski's geometry in a form 
similar to the standard differential geometry. 

Let us consider an auxiliary d-space (W, W), W = Gen(/3i, 02, ■■■,^5), where 

Pi-.W — s^R, f3i(t,x,y,z) =7TjO F(t,x,y,z), i = l,2,...,5. 
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This space is not diffeomorphic to the d-space with d-boundary (F(W),£p) w ^), since 
the function F: W — > R 5 is not one to one. Additionally, the generators Pi do not 
distinguish the following points p G dW := W — W°. Therefore, the d-space (W, W) 
is not a Hausdorff space. With the help of this space one can build a d-space with 
d-boundary (W,W) diffeomorphic to (F(W), 8f {w) ). 
Let qh be the following equivalence relation 

Vp, q : PQ H q & V/3 G W : (3(p) = (3(q). 

The quotient space W = W/pn can be equipped with a d-structure W := W/pn 
coinduced from W 

VV:= Gen(ftJ 2 ,...J 5 ), 

where 

A: W — > R, := ft(p), i = 1,2,..., 5. 

The symbol [p] denotes the equivalence class of a point p E W with respect to qh- 
The pair (W,VV) is a d-space [22 [231 [22]. 
Let us define the following mapping 

F: W — > F{W), F([p]) := F(p). 
Theorem 3.1 The d-space (W, W) is diffeomorphic to the d-space with d-boundary 

7(5) \ - 



Proof. The mapping F: W — > F(W) is a bijection. In addition, W is by con- 
struction a Hausdorff topological space with the topology given by the generators 
C211I31 ?]. The mapping 

F{\p\) = fa®, h(p),..., p s (p)) = {M\p\)M\p\),---M\p\)) 

is a diffeomorphism of the d-space (W, VV) onto its image (F(W), Sp}^\) = {F{W),£p^ w ^) 
[25]. □ 
According to Definition I2.3[ (W, VV) is the d-space with d-boundary for the flat 
FRW d-manifold (W°, W°), where the d-boundary is represented by the set djW := 

w-w°. 

The d-space with d-boundary (W, W) has been constructed with help of the d- 
space (W, W) and the relation qh- Generally speaking, d-spaces of the type of (M, C) 
and {M / qh,C / qh) have a lot of common features because of the isomorphism of the 
algebras C i C/qh- In particular, moduli of smooth vector fields X(M) and X(M/ qh) 
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are isomorphic [231 E2]- This property will, in the next parts of the paper, enable us 
to work with the help of the more convenient d-space (W, VV) instead of (W, VV). 

The d-structure W of the d-space (W, VV) is finitely generated by means of func- 
tions (3i, i = 1,2,. ..,5. However, in the next parts of the paper, we use a different, 
but equivalent, system of generators 

Oil := Pi -^2, "2 := Pi + (3 2 , Oii:=pi, i = 3, 4, 5. (3) 

Then the d-structure W has the form 

W = Gen(«i, a 2 , a 5 ). 

4 Time orientability 

The flat FRW model is a time orientable Lorentzian manifold M. By definition, there 
is a timelike directional field generated by a nowhere vanishing timelike vector field 
X. If X generates the directional field then the field AX generates it also, where A 
is a nowhere vanishing scalar field on M. The field X caries a part of information 
included in the casual structure of M, which enables us to define the direction of the 
stream of time and the succession of events [26] . 

The manifold structure and the casual structures of space-time are broken down 
at the initial singularity. In the hierarchy of space-time structures the Sikorski's 
d-structure is placed below the casual structure [15]. Therefore, the d-space with 
d-boundary (W, VV) of the flat FRW model is timeless independently of the fact that 
one of coordinates is called time and the moment t = is named the beginning of 
time. In this situation one cannot say that the cosmological singularity (d-boundary) 
is an initial or final state of the cosmic evolution. It is necessary to introduce a notion 
which would be a substitute of time orientability. 

Let (W°, W°) and (W, VV) be the pair of d-spaces described in the section [3j 
For convenience we can consider the d-space (W, VV) instead of the d-space with d- 
boundary (W, VV) according to the remark after Theorem 13. 11 In this representation 
the set of not Hausdorff separated points dW := W — W° is a counterpart of the 
d-boundary ddW° for the flat FRW d-manifold. 

Let in addition, X°: W° — > TW° be a timelike and smooth vector field without 
critical points, tangent to the manifold (W°, VV°), fixing the time-orientability on the 
manifold (W°,g). 

Definition 4.1 The d-space (W, VV) is said to be time oriented by means of a vector 
field X if 
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a) there is a nonzero vector field X:W — > TW tangent to (W, W) given by the 
formula 

I y\\Yii q ^ v X^{q){a\ w v) for p G dW, q <E 

b) and there is a function A G W, A(g) > (or A(g) < 0) for q G W° C W snc/i 
that the vector field V := AX is smooth on (W, W). 

v4 coordinate defined by means of X is called time and the moment t = t/je 
beginning of time t. We also say that the d-space (W, W) is oriented with respect to 
time t. 

In the flat FRW model ([1]) the cosmological time t is a time variable. The vector 
field of the form 

X °: W° — ► TW°, X°(p)(a°) := (4) 

where a G W°,p G W°, establishes the time orientation on (W°,g). The vector field 
is smooth on the manifold (W°, W°) since derivation X°: W° — > R w ° , X°(a°)(p) := 
X°{p)(a°) satisfies the condition X°(W°) C W° (Definition 



In the next parts of the paper, cosmological models for which the vector field X° 
can be extended on (W, W) are discussed. Then if remaining conditions of Definition 
14.11 are satisfied, (W, W) is a d-space with boundary dW of the flat FRW d-manifold 
which is time oriented with help of the following vector field 

X-.W-^TW, X(p)(a):=^M, (5) 

where a G W, p G W. In other words, the d-space (W, W) with boundary dW is 
time oriented with respect to the cosmological time t. 

Lemma 4.1 The mapping X:W — > TW , X(p)(a) := da(p)/dt is a vector field 
tangent to (W, W) if and only if for every t G D a , a(t) is finite and a(t) ^ 0. ■ 

Proof. A vector field is tangent to (W, W) if its value X(p)(a) is finite for every 
a G W and p G W. It is enough to check this property on the generators ai, «2, «5 
since the d-space (W, W) is finitely generated. Straightforward calculations show that 
the value of the field X(p)(a 2 ) is finite forpeW iff a{t) is finite and d(t) ^ for 
t G D a . Then, the value of X on the remaining generators is always finite. □ 
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Lemma 4.2 If V := XX, A G W, is a smooth vector field tangent to the d-space 
(W, W) and for every p G W° the value of the function A(p) 7^ then A(p) = /or 
p G cW. In ot/ier words, the smooth vector field V has a critical point at the boundary 

dw. m 

Proof. Every smooth function 7 G W is a local W-function on W (see Definition 
lA.lj) . This means that for every p e W there is C/ p G r w and / p G C°°(IR 5 ,IR) 
such that 7(g) = fp(cti(q), a 2 (q), ^(g)) for g G £/ p . Therefore, the value of 7 
in p = (0, x, y, z) G dW, where x, y, z are any, is a constant function of x, y, z: 
7(0, x, y, z) = fp(ax(p), a 2 (p), a 6 (p)) = / p (0, 0, 0, 0) = const. 

Now, let us suppose that X(p) 7^ for p G cW also. Smoothness conditions for 
the field V have the form: V^ai) G W,i = 1,2,..., 5. In particular, the conditions 
V(cti) G W and V(a 2 ) G W lead tod G W and a 7^ for p G W 7 . Then, for example, 
the function ^ := V(a 3 )/Xa, rj(p) = x, is a smooth function (77 G W) since A G W 
and A(g) 7^ for g G W 7 . This is a contradiction since this function is not a constant 
function on dW and therefore it is not generated by means of «j, i = 1, 2, 5. □ 



5 A smooth evolution with respect to cosmological 
time 

The vector field X defined by formula ([5]) is not smooth on (W, W). Its value on 
a smooth function is not necessarily smooth. This means, that there are functions 
a G W such that ^ W, where the derivation X: W — )■ IR 1 ^ is given by the 

formula X(a)(p) := X(p)(a). But also it can happen that, for a class of smooth 
functions, values of X on functions from this class are smooth. 

According to the definition of X, its restriction X° = X|n/o is a smooth vector 
field on (W°, W°), and therefore its value := X°(a°) = a on the smooth function 
a\ := ai| pi/0 is smooth: a \ G W°. The generator ct\ = a is by definition a smooth 
function on (W, W). According to the earlier argumentation, the value of X on the 
function a\, di := X(a\) = a, is not necessarily smooth in the Sikorski's sense. But, 
from the physical point of view, it is natural that the velocity of the expansion of the 
universe is a smooth function even at the beginning of time t = since a moment 
later, it is smooth ( 6t[ G W° ) both in the classical sense and Sikorski's sense. Let 
us distinguish a class of cosmological models with such a property. 
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Definition 5.1 An evolution of cosmological model ([T]) is said to be smooth from the 
beginning of cosmological time if 6l\ := X{a{) = a G W and a(t) ^ for t G D a . We 
shall also say that the cosmological model is smoothly evolving. 

Theorem 5.1 If the flat FRW model is smoothly evolving, then the d-space (W, W) 
with boundary dW of this model is time oriented with respect to the cosmological 
time t. A smooth vector field V defining time orientability on (W, W) has the form 
V:W-^W, V:=a\d/dt. ■ 

Proof. Proof consists of verification whether the following inclusion is satisfied 
V(a t ) = al^e W, z = l, 2,..., 5. □ 



6 The simplest smoothly evolving models 

Lemma 6.1 A function tp: W — > 1R which depends on the single coordinate t only 
is smooth (f G W if and only if, for every p G W , there is a neighbourhood U p G ryy 
and a function f G £^ such that <p\u p = f(ai,^)\u p , where 

£ : = axa2 - a \ - a\ - a 2 5 , £{t,x,y, z) = a(t) (6) 

./n a T 



Proof. Let a function ip be a smooth one. The d-structure W is finitely generated. 
Therefore for every p G W there are a neighbourhood U p G tw and a function 
y(5) g g(5) jjofifc ^ = J( 5 )(q/ 1; a2; Q/g)]^. Since is a constant function 

with respect to x, y, z, it has to be generated by means of such combinations of the 
generators admissible in the algebra W that a resulting function depends on t only. 
The evident form of the generators leads to conclusion that such combinations have 
the following form ct\ and £ := ct\Ct2 — ol\ — ol\ — ct\. Therefore /^ 5 ^ has the following 
structure /^(ai, a 2 , a 5 ) = a±a 2 — a\ — a\ — af) = f(a\, £), where / G 
The implication in the opposite direction is obvious. □ 
According to Lemma (16. 1ft . if di G W then, in a neighbourhood C/ p of p G W, 
the function di is of the form a>x\u p = f(&i,£)\u p - Especially interesting points in 
W are those not Hausdorff separated points from the boundary p, G dW, p, = 
(0,x,y,z). For a given /, from the whole class of all possible neighbourhoods U Pm 
such that di|(7 p> = f(oti,£)\u p . one can select the maximal neighbourhood in the 
sense of inclusion. This means in practice that the maximal neighbourhood has the 
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form Up, = D a x M 3 , where the domain D a of the scale factor is a subset of the set 
{t G [0, oo) : ^ a(t) < oo, < a(t) < oo} such that G D a . In the neighbourhood 
Up, the smoothness condition for a\ (see Definition I5.ip has the form 

a(t) = f{a(t),a(t) / (7) 
Jo a ( T ) 

where 

/(0,0)>0, /g£ (2) . (8) 

Formula (JTj) is an equation for a(t) with an initial condition a(0) = 0. The function 
/, is in a principle, arbitrary. The only restriction on / is condition (JBJ) which is a 
consequence of Lemmas 14.11 16.11 and the physical assumption that the real universe 
expands from the initial singularity. The simplest choice is the following function 

f(zi,z 2 ) :=p + 7i-2i + 72^2, 0>O, 7i,72 eR- 
Now, the smoothness equation (JTj) for a(t) has the form 

a(t) = P + lia{t)+ l2 a{t) / — -. (9) 

Solutions of ((9]) have to satisfy the following conditions 

o(0) = 0, d(0) =P>0, a(t) > 0, a(t) > for t > 0. (10) 

Proposition 6.1 When 72 = t/ien solutions of the smoothness equation (Tj|) satis- 
/ym^ conditions ICTty) have the form 



a(t) = (3t, te [0,oo), /or 7l = 0, (11) 

and 

a(t) = ^(e^ - 1), t G [0, 00), /or 7l ^ 0. (12) 
7i 



Proof. Obvious calculus. □ 
Solution (ITT]) represents the well known model of the universe which expands 
with the constant velocity a — f3 and which is a solution of the Friedman's equations 
with the following equation of state p = —g/3. For 71 > solution ( Tl2]) is the 
universe model which is asymptotically ( t — > 00) the de-Sitter model. The model 
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expands from the very beginning with a positive acceleration. The parameter 71 can 
be asymptotically interpreted as a cosmological constant. When 71 < cosmological 
model ffl~2]) describes an expanding universe, and the expansion slows down from the 
very beginning. For great t, the size of universe fixes on the level a(t) ~ lim^oo a(t) = 
0/\"fi\ and a and a tend to zero when t — > 00. Such a universe asymptotically becomes 
the Minkowski space-time. 

In the case 72 7^ 0, let us introduce the following auxiliary symbols 

7i : =7i72/V3, 72 := sgn(7 2 )/372/3, if := 72 > 0, 

a(if) := ^® , *(if ) := 72 *(*0/V3, (13) 
where H(t) := a(t)/a(t). 

Proposition 6.2 1/72 > inen solutions of smoothness equation (j5J) have the form 

a(K) = (K-j 1 - arccothK)- 1 , t(K) = l°° a ^ y J ydy } (14) 

Jk r - 1 

where if G (iff, 00) , and if/ is a solution of the following equation 

Kf — arccothifj = ji, if/G(l,oo). (15) 

■ 

Proof. Solution of an elementary differential equation. □ 
Proposition 6.3 7/ 72 > then 

1- if li ^ 0, acceleration d(t) > for t G (0, 00) , 

2. if ji < 0, acceleration a(t) < for t G (0, i*), a(i*) = and d(t) > /or 
i G (£*, 00), where t* := t(if*) and if* zs a solution of the following equation 

if 

71 + ^ 2 1 + arccothif* = 0, if* G (if/, 00). (16) 



Proof. By obvious calculation. □ 
In the case 72 > 0, there are two essentially different scenarios of a smooth evolu- 
tion with respect to the cosmological time t: 
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a) If 7x ^ the model accelerates from the very beginning and expands indefinitely. 

b) In the case 71 < 0, initially the expansion slows down, but at the moment 
the unlimited and infinitely long accelerated expansion is initiated (see Figure 



a(t) for 72 > 




Figure 1: Scale factor a(t) for smooth solutions with 72 > 0. When 71 ^ models 
expand with acceleration. For 71 < solutions initially expand with a negative accel- 
eration but at a moment, denoted by a circle on the graph, an accelerated expansion 
is initiated. The black point denotes the initial singularity. 



Proposition 6.4 If j 2 < 0, solutions of smoothness equation Qj have the following 
form 

a(K) = (ir-^i-arctanK + vr^)- 1 , t(K) = f°° (17) 

Jk y + 1 

where K 6 (Kf, 00) and Kf is a solution of the following equation 

# / -arctanif / + 71-/2 = 71, K f e R. (18) 

Proof. Solution of an elementary differential equation. □ 

Proposition 6.5 If 72 < and 71 < tt/2, then cosmological model ([I]) has a final 
curvature singularity at t s < 00. The set [0,t s ) is a domain of the scale factor a(t), 
where t s := t(0) and the function t(K) is given by formulae ( fT7|) and ( Tl3l) . 
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Proof. Proposition is the result of a fact that some of components of the curvature 
tensor are undefinite at t s , because a(t) — > — oo when i — > i~. □ 
According to Propositions 16.21 and I6.4[ solutions of smoothness equation fl9]) are 
defined in the domain (Kf, oo). But in the case of solutions discussed in the Propo- 
sition [631 there appears an additional restriction for the domains of a(K) and t(K) 
coming from the geometrical interpretation of a(K) as the scale factor for a flat FRW 
model. The final curvature singularity ends the evolution of the model. 

Proposition 6.6 1/72 < then 

1- if lx ^ 7r /2 ; acceleration d(t) > for t G (0, oo), 

2. if < 71 < tt/2, the scale factor a(t) has the inflexion point at i* (a(t%) = 0), 
a{t) > fort E [0, £*) and d(t) < for t G (t*, t s ), 

3- if lx ^ 0, acceleration d(t) < for t G [0, t s ), 

where t* := t(K*). The quantity K* is a solution of the following equation 

— — — - — arctanK* + 7r/2 = 7!. (19) 
Ki + \ 



Proof. By properties of the function given by formulae (I17p . □ 
If 72 < 0, the smooth evolution of the universe with respect to the cosmological 
time t can proceed on three different ways: 

a) If 71 ^ 7r/2, a smooth accelerated evolution starts from the initial singularity. 
The acceleration goes on continuously during an infinite period of time. 

b) For 71 G (0,7r/2), these smoothly evolving universes initially accelerate but 
the acceleration is slowing down so as to change, at t — i*, into a deceleration. 
Smooth evolution ends at the final curvature singularity within the finite period 
of time [0, i s ). These models have two singularities: the initial and final one. 

c) Models with 71 ^ start their evolution in the Big-Bang and decelerate. The 
rate of expansion is strongly slowing down and these models end their evolution 
at curvature singularities in a finite time t s . Models of this class have also two 
singularities. 
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a(t) for 72 < 



a 



0.5 



2.5 



2.0 



0.0 



1.5 



1.0 




0.0 



t 



Figure 2: Scale factors a(t) for smoothly evolving models with 72 < 0. Curves on the 
graph with 71 ^ tt/2 represent accelerated solutions. For models with 71 G (0,7r/2) 
initially accelerating expansion is slowing down. At the moment (circles on the 
graph) the acceleration changes into a strong deceleration. When 71 ^ solutions 
expand with negative acceleration. The black points on the graph denote initial and 
final curvature singularities. 

7 Primordial matter and dark energy 

Let us assume that solutions of the smoothness equation represent cosmological mod- 
els. Then the Friedman equations 



can serve as a definition of a pressure p and energy density g of a kind of cosmological 
fluid which causes the smooth evolution of models, where p = up and g = Kg. In the 
present paper this fluid is called the cosmological primordial fluid, or primordial fluid 
for brevity. 

Proposition 7.1 If '71 G K and 72 = then the equation of state for the primordial 
fluid has the form 



p = —2a/a — a 2 /a 2 , 



g = 3d 2 /a 2 , 



(20) 




(21) 



In addition 
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a) the energy density g is a decreasing function of cosmological time and lim^o+ Q(t) 
oo, 

b) at the initial moment p(0) = \im t ^ + p(t) = — oo, and the remaining details of 
the dependence p(t) are shown in Figure® 



c) 



72 = 


Q 


V 


SEC 


WEC 


7i > 


lim^oo g(t) = 3j[ 


lim^oo p(t) = -37i 


no 


yes 


7i ^ 


lim^oo g(t) = 


]xm t -+<x>p(t) = 


yes 


yes 



Proof. An elementary calculus. □ 
The abbreviations SEC and WEC on the above and next tables denote the strong 
and week energy conditions and the statements below are answers to the question of 
whether the strong or the week energy conditions are satisfied. 



p(t) for 72 = 




Figure 3: Dependence pit) for 72 = 0. In the case 71 < the pressure has a single 
positive maximum. 



In the case 72 7^ it is convenient to introduce the following abbreviations 



P-=P% 2 , Q-=Ql2 2 - 



(22) 
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Proposition 7.2 If ^2 > and 71 G R £/ien the equation of state of the primordial 
fluid has the form 

2 1 2 

p = -- - -g - -(7! + arccothv/^)(A/^ - 1/ a/^). (23) 

In addition 

a) energy density g is an increasing function of time and lim t ^ + Q{t) = +00, 
lim^oo g(t) = Kj where Kf is a solution of equation (|15p . 

b) pressure at the beginning and end of the evolution is lim t ^ + = — oo ; \im t ^. oc p(t) = 
—Kf, and the remaining details of dependence p(t) are shown in Figure® 

c) the week energy condition is satisfied during the whole evolution, 

d) for 71 ^ the strong energy condition is broken down, 

e ) if li < 0, the strong energy condition is satisfied for t e (0,t*). For remaining 
t > t* this condition is broken down. The moment t* is defined in Proposition 

Proof. An elementary calculus. □ 

p(i) for 72 > 







71 < 






i 


71 =0 






71 ^ . 







1 2 3 7 4 



Figure 4: For ji > —1.35, pressure is a decreasing function of cosmological time. For 
remaining 71, function p(t) has both the maximum and minimum. The minimum is 
not well visible on the graph. 
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Proposition 7.3 If 72 < and 71 G K i/ien £/ie equation of state of the primordial 
fluid has the following form 

2 1 2 r r j— 

P = +3 - gff- 3(71 - arccotV^)(V^+ 1/Vq). (24) 

energy density g is a decreasing function of the cosmological time, and lim i _ >0 + = 
+oo ; 

6^ initial pressure is p(0) := lim t ^ +p(^) = —00, and t/ie remaining details of 
dependence p(t) are shown in the Figure El 



c) 



72 < 




p 


SEC 


WEC 


7i ^ tt/2 


lim^oo g(t) = K'j 


lim^ 00 p(t) = -K'j 


no 


yes 


< 71 < tt/2 


\im t ^ ts g(t) = 


]im t ^ ts p(t) = +00 


no/yes 


yes 


7i ^0 


limt-rt, = 


lim^-rt, p(t) = +00 


yes 


yes 



d) if < 71 < 7r/2, t/ie cosmological fluid violates the strong energy condition for 
t G (0,t*). In i/ie remaining range t G (t*,t s ), the SEC zs satisfied. 

Quantities Kf, t s and t* are defined in Propositions 16.4116.51 and 16.61 

Proof. Elementary calculations. □ 



The simplest solution, 71 = 72 = 0, of smoothness equation (|9]) represents a 
model filled with the primordial fluid with the equation of state p = —l/3g. In the 
present paper, this fluid is called a 70-matter. In the case of the following parameters 
system {7! 7^ 0,7 2 = 0}, the primordial fluid consists of the 70-matter enriched by a 
material ingredient connected with the generator a\(t) = a(t) in formula fl9]). This 
enriched primordial fluid we call a 71-matter when 71 < 0, or a 71-energy when 
71 > 0. Similarly, in the case of the following parameter system {7! = 0, 72 7^ 0}, the 
primordial fluid composed of the 70-matter and a matter connected with the generator 
Q2, through the function £ in formula fl9]) , we call a 72-matter when 72 < 0, or a 
72-energy when 72 > 0. 

Taking into account Proposition 17. II 71 -energy has the properties of a dark energy. 
This energy causes the expansion to accelerate. During the evolution acceleration 
grows to infinity. After an infinitely long evolution pressure and energy density reach 
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71 ^0 



p(t) for 72 < 



10 



P 



< 7x < tt/2 




5 



71 tt/2 



-5 



3 



4 



5 



t 



6 



Figure 5: The graph shows the great qualitative differences in the p(i) dependence 
for various ranges of 71. 

the finite values Pf = —Sjf and p/ = 37^ respectively. Let us notice that then the 
following equation of state for the cosmo logical constant is satisfied, pf = —Qf- The 
71-energy reaches this property at last stage of the evolution. 

In contrast to 71-energy, 7 1 -matter satisfies the strong energy condition during 
the whole evolution. It has an attraction property. Therefore the expansion of the 
cosmological model is slowing down in such a manner that at the last stages of the 
evolution the model becomes static. The 71-matter changes its properties during the 
evolution. Initially, it has a negative pressure. But later it transforms itself into 
a kind of matter with a positive pressure. At the last stages of the evolution the 
pressure and the energy density of 71-matter become zero: pf = and Qf = 0. After 
an infinite period of time since the Big-Bang this cosmological model becomes, in an 
asymptotic sense, the Minkowski space-time. 

A model of the universe filled with 72-energy monotonically accelerates. Initially, 
the jostling property of 72-energy has a small influence on the expansion but its 
inflationary power is disclosed at the last stages of evolution of the model. The 
initially infinite energy density strongly decreases and at the end of the evolution 
is on the level of Qf = ^\Klp where K f w 1.19 is a solution of equation f lT5|) for 
71 = 0. A negative pressure rapidly grows from —00 to a finite level of pf = —j^^of- 
At the end of the evolution the equation of state is as for the cosmological constant: 
Pf = —Qf- During the whole evolution the strong energy condition is violated. 72- 
energy can be interpreted as a dark energy of different kind then 71-energy. Details 
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can be found in Figures [TJ @] and in Proposition 17.21 

72-matter has a strong attraction property and therefore the expansion of a model 
with such a fluid is rapidly slowing down. Acceleration quickly decreases from to 
— oo in a finite period. At the and of the evolution the model is stopped a(tf) = 0, and 
its scale factor reaches the maximal, finite value. The final curvature singularity ends 
the evolution. Properties of 72-matter are changing during the evolution. Pressure 
rapidly increases from —00 to +00 in the finite period [0, tf). Simultaneously, the 
energy density decreases from +00 to zero independently of the fact that the scale 
factor is finite (a(tf) < 00) at tf. This kind of matter has very interesting properties 
at the end of the evolution: it has a slight energy density but simultaneously a huge 
positive pressure. Details can be found in Figures HJ and in Proposition 17.31 

A mixture of ^-energy and ^-energy (7! > 0, 72 > 0). During the whole period 
of evolution of this model the mixture has the properties of a dark energy. Both 
components interact with each other causing increased acceleration. In the last stages 
of the infinitely long evolution, the equation of state for the mixture has the form 
of equation of state for the cosmo logical constant f?/ = ~Pf = 72-^/5 where Kf is a 
solution of equation ( TT5|) . 

A mixture of 'ji-matter and 72-matter (7! ^ 0, 72 < 0). This kind of the primordial 
fluid satisfies the strong energy condition during the whole finite period of evolution. 
The mixture is a fluid with interacting components. At the final singularity 72-matter 
absorbs, in a sense, 71-matter and finally the mixture vanishes, Qf = at an infinite 
pressure. An admixture of 71-matter into 72-matter shortens the life time of the 
cosmological model. 

A mixture of •yi-matter and ^-energy (71 ^ 0, 72 > 0). Components of the 
mixture are interacting fluids. The beginning of the model evolution is dominated 
by 71-matter. The universe expands with a negative acceleration and the mixture 
satisfies the strong energy condition. But the influence of 72-energy is still rising. 
For t > the evolution is dominated by 72-energy. Since t = t*, the mixture 
has properties of a dark energy and changes the further evolution into accelerated 
expansion. During final stages of the evolution the equation of state has the following 
form Qf = —pf = 72-^/) where Kf is a solution of equation (fT5|) . 

The evolutionary behaviour of the model is extremely interesting because such 
an evolution is qualitatively consistent with the observational data of la type super- 
novae PQ[2]. Preliminary quantitative investigations of the consistency of the discussed 
model have been carried out with the help of the H Q b s (z) dependence published in [27]. 
Results of the best-fit procedure depend on H . For 70.6 ^ H ^ 77.8 kms _1 Mpc _1 
[28] the best-fit parameters are in the range 71 G [-1.829, -1.075], 72 G [2.82, 4.011] 
xl0 _4 Mpc _1 and Xmin e [8.66, 9.76] (see also Figure^]). Values of parameters 71 and 
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72 enable us to find: an age of universe to £ [13.561, 14.241] x 10 9 y, the moment of the 
acceleration beginning t* G [7.427, 8.471] x 10 9 y, the Hubble constant in the acceler- 
ation moment H* G [106.31, 108.86] x kms^Mpc" 1 and redshift z* G [0.648,0.743], 
where 

z(K) := ^ - 1, t Q = c-H{K ), U = c-H{K*), K := ^H . 
V3(3a(K) C7 2 

Quantities K, a(K), t(K) are given by formulas (lT3i) and (Til]) , and if* is a solution 
of equation ( [TBI . In the present Section bar over quantities denotes that we use the 
systems of units in which c ^ 1. 




0L ■ : z * - . ■ J 

0.0 0.5 1.0 1.5 2.0 2.5 

z 

Figure 6: Comparison of the theoretical H(z) dependence for the smoothly evolved 
model (71 < and 72 > 0) with H b s (z) for H = 72.0 km s _1 Mpc _1 . The solid 
line represents the smoothly evolved model with the best-fit parameter values 71 = 
-1.241, 72 = 0.000312MPC" 1 and xLn = 8-79. In this case t = 14.103 x 10 9 y, 
U = 8.241 x 10 9 y, H* = 106.87 kms^Mpc- 1 and z* = 0.668. The dashed line 
represents the best-fit of the ACDM model. 



A mixture of 71- energy and 72-matter (71 > 0, 72 < 0). For t G [0,£*) the mixture 
has properties of a dark energy and therefore this model accelerates from the very 
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beginning. But later, 7 2 -matter component begins to play a bigger and bigger role. 
The further evolution depends on the value of the parameter 71. 

When 71 e (0, 7r/2), the repulsive properties of 71-energy are not able to dom- 
inate the evolution and the acceleration is gradually stopped because of a greater 
and greater attractive influence of 72-matter. The moment t*, defined in the Propo- 
sition 16.6} is the end of acceleration. Starting from this moment, the expansion is 
strongly slowing down till the final singularity. The behaviour of pressure is interest- 
ing, Figure El for t > i*. After t* the mixture behaves as a dust (p«0). But later, 
pressure rapidly increases to infinity at the final singularity. Simultaneously, because 
of expansion, the energy density decreases to zero. 

When 71 ^ 7r/2, repulsive influence of 71-energy is dominating during the whole 
period the infinitely long evolution. The mixture acts like a dark energy causing 
acceleration, independently of the 72-matter presence. As for previously considered 
accelerating models, the equation of state, in the last stages of the evolution, is as 
the one for the cosmological constant, i.e., Qf = —pf = l\K 2 ^ where Kf is a solution 
of equation f|T8|) . 

8 Smoothly evolved models in a neighbourhood of 
singularity 

The function / appearing in the smoothness equation (J7J) can be expanded into a series 
in a neighbourhood of the point (0, 0): f(zi, z 2 ) = (3 + d\f(0, 0) ■ Zi + d 2 f(0, 0) ■ z 2 + .... 
Then the smoothness equation assumes the following form 

/"* dr 

a(t) = (3 + d 1 f(0, 0) • a{t) + d 2 f(0, 0) ■ a(t) / — - + ... . 

Jo a(r) 

For sufficiently small t, or equivalently in a small neighbourhood of the initial singular- 
ity, one can omit higher powers of the expansion and consider the smoothness equation 
in the linear approximation. If one assumes that 71 := <9i/(0, 0) and 72 := (92/(0,0), 
the above equation, in the linear approximation, is identical with the smoothness 
equation OH]) for the simplest smoothly evolved models (Section EJ). 

The above observation leads to the conclusion that properties of the solutions 
studied in Sections |6] and for small t, are typical for every smoothly evolving 
flat FRW model in a small neighbourhood of the initial singularity. In particular, 
every smoothly evolving model during the initial stages of its evolution is filled with 
a cosmological fluid which is 7j-matter or 7 r energy, % = 1,2, or one of mixtures 
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described in Section [7J These fields satisfy the following approximate equation of 
state 

Q 2 7l _ 4_ 2 2 7l 7 2 2 1 2 4 7| 3 4 7 2 6 
3V3 3 V3 y/Q 9 45 g 2 

where 71 G R, 72 ^ 0, e = sign(7 2 ). 

9 Summary and Discussion 

a) The discussion of the equivalence principle and its breaking down has led to the 
formulation of the d-boundary notion which is, roughly speaking, an "optimal" set 
of points at which the equivalence principle is broken. A space-time d-manifold with 
the attached d-boundary is not only a topological space but also an object with a 
rich geometrical structure called differential space (see Appendix |AJ . The effective 
construction of both d-boundary and d-space with d-boundary for any space-time 
d-manifold has also been described (see Section [2]). 

b) For every flat FRW cosmological model with the initial singularity the d-space 
with d-boundary (W, VV) has been constructed. In this case the d-boundary is a 
single point (see Section [3]). 

c) In the d-space formalism it is possible to extend the concept of time orientability 
to the d-boundary of the FRW models. In this way, the intuitive understanding of the 
beginning of the cosmological time obtains a precise mathematical form (see Section 
H|). However, not every flat FRW model with the initial singularity has the well 
defined beginning of the cosmological time. 

d) In the whole class of flat FRW models with the well defined beginning of 
the cosmological time we distinguish the large subclass of the so-called smoothly 
evolved models (see Definition 15. ip . The condition defining this subclass is called the 
smoothness equation. The most practical form of this equation is given by formula 

0- 

e) The simplest flat FRW models with the well defined beginning of the cosmo- 
logical time have been found in the explicit form. This set of models can be divided 
into six qualitatively different classes (Lemmas 16. 1[ 16721 and \QA\ i. The most important 
classes are 

• solutions with parameters 71 < 7r/2 and 72 < which have two curvature sin- 
gularities: the initial singularity and the final singularity (see Figure |2]), 

•• the subfamily with 7l < and 72 > 0, being qualitatively consistent with the 
observational H b s (z) data from [27]. The quantitative consistency with the data 
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is discussed in Section [7] and is shown in Figure EJ The level of the consistency 
is similar to that for the ACDM model. 

f) The Friedman equation without cosmological constant, in application to the sim- 
plest solutions, may serve as a definition of energy density g and pressure p of a 
cosmological fluid (primordial fluid) which causes the smooth evolution. This strat- 
egy enables us to find main phenomenological properties (in particular, the equation 
of state) of this primordial fluid (see Lemmas 17. 1[ 17.21 i 17. 3p . 

g) In Section [7] we present an interpretation of the primordial fluid as a mixture 
of more elementary interacting primordial fluids 7j-matter and 7j-energy, % = 1, 2. 

h) From the analysis of smoothness equation (J7J) in a neighbourhood of the initial 
singularity, i.e for small t, one can conclude (see Section [8]) that in earlier stages of 
the evolution every smoothly evolving flat FRW model is filled with a primordial 
cosmological fluid with properties characteristic for particular solutions found in this 
paper (Sections Eland [7]). 

i) It is very surprising that without any assumption concerning physical nature of 
the cosmological fluid , among the simplest solutions of the smoothness equation (J7J), 
it is possible to find models consistent with the observed evolution of the Universe 
(see Figure [6] ) . The existence of the well defined beginning of the cosmological time 
was the only requirement which has been assumed. Primordial mixture of fluids 71- 
matter and 72-energy is a consequence of this simple assumption. However, we cannot 
expect that every detail of the cosmic evolution can be determined in this way. The 
material ingredients such as radiation, dust and dark matter should also be taken 
into account. 
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A Sikorski's differential spaces 

Let C be a non empty family of real functions defined on a set M. The family C 
generates on M a topology denoted by the symbol tq. It is the weakest topology on 
M in which every function from C is continuous. 

Let A C M be a subset of M and let the symbol C\A denotes the set of all 
functions belonging to C restricted to A. On A one can define the induced topology 
r c H A = t c \a- The topological space (A, t c \a) is a topological subspace of (M, r c ). 
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Next, we introduce two key notions in the d-spaces theory: a) the closure of C 
with respect to localization and b) the closure of C with respect to superposition with 
smooth functions from S (m) := C°°(M m , R), m = 0,1, 2, .... 

Definition A.l A function 7: A — > R is said to be a local C -function on a subset 
A C M if, for every p G A, there is a neighbourhood U p G tq\a and a function <fi G C 
such that j\U p = <f>\U p . The set of all local C-functions on A is denoted by Ca- 

It is easily to check that in general C\a C Ca- In particular C C Cm- 

Definition A. 2 A family C of real functions on a set M is said to be closed with 
respect to localization if C = Cm- 

Definition A. 3 A family of functions C is closed with respect to superposition with 
smooth functions from £^ m \ m = 0, 1,2, if for every function uo G and for 
every system of m functions <pi, <f2, ■-, <fm G C, the superposition ou(ipi, <f2, <fm) is 
a function from C; ou(<pi, (fi2, <p m ) £ C. 

The above described system of concepts make possible to define an object (a d-space) 
which is a commutative generalisation of the d-manifold concept. 

Definition A. 4 A pair (M,C), where M is a set of points and C a family of real 
functions on M , is said to be a differential space (d-space for brevity) if 

1. C is closed with respect to localization, C = Cm, 

2. C is closed with respect to superposition with smooth functions from £( m ) . 

The family C is called a differential structure on M (d-structure for brevity) and the 
set M a support of the d-structure C. Functions if G C are called smooth functions. 

Every d-space is simultaneously a topological space with the topology r c given 
by the d-structure C in the standard way. The d-structure itself, with the usual 
multiplication, is a commutative algebra. The notion of smoothness, given by the 
condition ip G C, is an abstract generalization of the smoothness notion for functions 
defined on R™. Differential structure, by definition, is a set of all smooth functions 
on M. There are no other smooth functions on M. This class of smooth functions 
may consists of functions which are not smooth in the traditional sense. This is a 
great advantage of the d-spaces theory. The simplest example of a d-space is the 
n-dimensional Euclidean d-space (W,S {n) ), where S {n) = C°°(W n ,R). 
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There exists a procedure to construct a d-structure with the help of a chosen set 
of real functions on M. Let us denote it as Cq. The method consists in adding to 
a given Cq missing functions so as to satisfy the axioms of the closure with respect 
superposition with smooth functions and the closure with respect to localization. The 
closure with respect to superposition with smooth real functions on M n , is denoted 
by mathematicians by sc(C ) and the closure with respect to localization is denoted 
by (Cq)m (see Definition lA.lj) or lc(C ). 

It is easy to check that 

Lemma A.l Let Cq be a set of real functions defined on a set M. The family of 
functions C := lc(sc(Co)) = (sc(C ))m is the smallest, in the sense of inclusion, d- 
structure on M containing Cq. 

Sometimes one uses the following abbreviation: C = Gen(C ) := lc(sc(C )) = (sc(C ))m- 

Definition A. 5 The set Cq in Lemma \A.l\ is said to be a set of generators. Functions 
V? G Cq are called generators of the d-structure C := lc(sc(Co)). If Cq is finite then the 
d-structure C is called finitely generated. 

The method of constructing a d-structure with the help a set of generators is the 
great advantage of the Sikorski's theory, especially in the case of finitely generated 
d-spaces such as, for example, the d-space with singularity associated with the flat 
FRW world model. 

Definition A. 6 If (M,C) is a d-space and A C M then the d-space (A,Ca) is said 
to be a differential subspace (d-subspace) of the d-space (M, C). 

The above definition enables us to determine a d-structure for any subset A of M. 
It is enough to "localize" every function from the d-structure C to A. In the case of 
a finitely generated d-spaces (M, C), a simpler situation occurs. Then C = Gen(C ), 
Cq := {/?i,/3 2 , ■ ■■,0 n } , n 6 N, where 01,02, ...,(3 n are given functions. The d-structure 
Ca is given in terms of generators Cq = Cq\A which is a set of restrictions of the set 
C to A. Then C A = Gen(C ). 

Definition A. 7 Let (M,C) i (N,V) be a d-spaces. 

1. A mapping f: M — > N is said to be smooth if\/0&T): o f 6 C. 

2. A mapping f: M — > N is said to be a diffeomorphism from a d-space (M, C) to 
a d-space (N, T>), if it is a bisection from M to N and both mappings f: M — > N 
and f~ l : N — > M are smooth. 

In this case, we say that (M, C) and (N, T>) are diffeomorphic. 
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A smooth mapping / transforms smooth functions on N onto smooth functions 
on M. The notion of d-spaces diffeomorphism is the key notion from the point of 
view of the present paper. If there is a diffeomorphism / between d-spaces (M, C) and 
(N,T>) then these d-spaces, from the viewpoint the d-spaces theory, are equivalent. 

Definition A. 8 Let (M,C) i (N,V) be d-spaces. A d-space (M,C) is said to be 
locally diffeomorphic to the d-space (N, D) if for every p G M there is U p G tq and 
a mapping f p : U p — > f p {U p ) G t v such that f p is a diffeomorphism between the d- 
subspaces (U p ,Cu p ) and (f P (U P ),'^ > fp(u p )) of the d-spaces (M,C) i (N,T>), respectively. 

Definition A. 9 A d-space (M,C) is said to be an n-dimensional d-manifold, if it is 
locally diffeomorphic to the d-space (M. n ,£^). 

Applying Definition IA. 81 to Definition \A.9\ leads to condition (*) in Definition 12.11 
Local diffeomorphisms f p are obviously maps. A set of maps forms an atlas. It turns 
out that Definition IA.9I is equivalent to the classical definition of d-manifold [T2] . 

Theorem A.l Let (M,C) be a finitelly generated d-space with the structure C gen- 
erated by a finite set of functions: C := (3 2 , •••> Pn}, C = Gen(C ). Then the 
mapping F:M — > M. n , F{p) := ((3i(p), (3 2 (p), (3 n (p)) is a diffeomorphism of the 
d-space (M,C) onto the d-subspace (F(M),£p) M .) of the d-space (M n ,£( n) ). 

Proof, see [25] . 

The d-subspace (F(M),£$ M) ) of (R n ,S^) is an image of the d-space (M,C) in 
the mapping F. Theorem IA.1I is called the theorem on a diffeomorphism onto the 
image and it is important for the construction of the d-space for the flat FRW model 
with the initial singularity. 

Definition A. 10 A mapping v: C — > R is said to be a tangent vector to a d-space 
(M, C) at a point p G M if 

1. Va, f3 G CVa, b G R : v(aa + b/3) = av(a) + bv(/3), 

2. Va, $ G C : v{ap) = v{a)(3{p) + a{p)v{(3), peM. 

The set of all tangent vectors to (M, C) at p G M is said to be a tangent vector space 
to (M,C) at p G M and is denoted by T p M. The symbol TM denotes the following 
disjoint sum: 
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Definition A. 11 Let (M, C) be a d-space. The mapping X:M — > TM such that 
\/p G M, X(p) G T p M is said to be a vector field tangent to (M,C). 

With help of a vector field X: M — > TM one can define the following mapping 

X:C — >R M , X(a)(p) := X(p)(a), 

where a G C i p G M. The mapping is linear and satisfies the Leibnitz rule Va, /3 G 
C : X(a(3) = X(a)(3 + aX((3). Therefore, it is a derivation and a global alternative 
for the definition of a vector field (1A.11I) . 

Definition A. 12 A vector field tangent to (M,C) is said to be smooth if the mapping 
X satisfies the condition: X(C) C C. 

Definition A. 13 Let f: M — > N be a smooth mapping. The mapping f* p : T p M — > T 3 
given by the formula 

Wv ET p M,f3eV: U(v)(f3) := v((3 o /), 

is said to be differential of the mapping f at the point p G M. 

Let us define the following mapping id a'- A — > M, idA{p) = p, where A C M. 

Definition A. 14 A vector field Y: M — > TM on (M,C) is said to be tangent to a 
d-subspace (A,Ca), if there is a vector field X: A — > TA on (A,Ca) such that 

VpeA: Y(p) = id A * P (X). 
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